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ABSTRACT 

We study a membrane - anti-membrane system in Matrix theory. It in fact ex- 
hibits the tachyon instabihty. By suitably representing this configuration, we 
obtain a (2+l)-dimensional U{2) gauge theory with a 't Hooft's twisted bound- 
ary condition. We identify the tachyon field with a certain off-diagonal element 
of the gauge fields in this model. Taking into account the boundary conditions 
carefully, we can find vortex solutions which saturate the Bogomol'nyi-type 
bound and manifest the tachyon condensation. We show that they can be 
interpreted as gravitons in Matrix theory. 
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1 Introduction 

It is well-known that a tachyon develops when a D-brane and an anti-D-brane come up to 
some critical interval of order of string scale |jl|, 0, |^ • Hence the system of a brane and an 
anti-brane is unstable, and the tachyon fields roll down to a stable point p. Recently Sen 
argued that the D-brane - anti-D-brane pair is annihilated via the tachyon condensation 
and in this process certain stable D-branes are produced as soliton solutions on this system 
0; i) 0! HI- These stable D-branes are basically classified into two classes. One is stable 
non-BPS D-branes with codimension one [^ |^, 0, ||, P]. They are stabilized by certain 
orbifold or orientifold projections, and can be identified with kink solutions. The other 
is BPS D-branes with codimension two, and they can be interpreted as vortex solutions, 
or kink solutions on unstable D-branes with codimension one, on the brane - anti-brane 
pair P, § . Both classes of stable D-branes has been understood in terms of the K-theory 



in more general way |Tl[, and the author in |0 suggested possible applications of this 
line of arguments to Matrix theory. 

In this paper we study the latter mechanism in Matrix theory |l^. We consider 



a specific example, a system of membrane - anti-membrane pair,0 in which gravitons 
are supposed to be produced as vortex solutions on the membrane worldvolume via the 
tachyon condensation: 



M2 + M2 ' '^ ; graviton (vortex), (1) 

where (T) denotes the vacuum expectation value of the tachyon fields. 
This corresponds to the mechanism in type IIA theory 0, 

D2 + D2 ^^^""^ > DO (vortex) or DO (anti- vortex). (2) 

In Matrix theory we should not expect the appearance of anti-gravitons (corresponding 
to DO in type IIA theory), for the matrix description of M-theory is thought of as an 
infinite momentum limit or as a light-like limit of M-theory. In fact we will see below that 
we cannot have an anti-vortex solution in the matrix calculation. 



2 The setup 

In string theory the tachyon potential can in principle be computed, but it is hardly 
possible to carry it out in all order. In Matrix theory, however, we in fact have a simple 
tachyon potential of quartic order, due to having taken a special limit of M-theory |l^ . 



■^Hereafter we will denote membranes and anti-membranes as M2 and M2 respectively. We will also 
use the abbreviation Dp and Dp respectively for Dp-branes and anti-Dp-branes. 



Now let us see how tachyon modes appear in the background of the M2-M2 pair in Matrix 



theory, which was studied in detail in |15, 16 1. The matrix theory action is 

^ = Wm^I"^^^' {^(AXT + l[X\Xf + vi/^A^ + ^^W[X., ^]} , (3) 

where X* (i = 1, ■ ■ ■ , 9) and \E' are NxN hermitian matrices. R is the radius of a light-like 
circle and Mp is the eleven dimensional Planck mass. We use a convention dt — i[Ao, ■] for 
the covariant derivative Df. 



We can express the M2-M2 pair spreading over the (8,9)-directions of interval b in the 
7-direction by the background matrices 



B' 



Pi 
P2 



fiS 



Qi 







-Q2 



B' 




b 



(4) 



where Pi and Qi satisfy the commutation relation [Pi, Qi] = —ic with a constant c which is 
a charge density of the membrane |]T^. The upper diagonal block corresponds to M2 and 



the lower to M2. The background matrices for the other directions, the gauge potential 
Aq and the fermions are vanishing. For later use we shall represent Pi and Qi by 



P 



-id^, + A^, 



Xil 1 



^Xi ^yi Oy. Ax^ 



Qi = c{-idy^+Ay^), (5) 

^ ^11^ . In this representation the trace in 



with the field strength F^^y. 

the action is replaced with integrals over Xi and ?/j. For our purpose, we are only interested 

in the type of a trace 

2 



Tr 



Ai X 

Xt D. 



Dy-. Y 

Y^ a, 



(6) 

where D^^ and Dy^ denote, in the representation given in eq. (^), covariant derivatives in 
which the gauge fields depend only on xi and yi for i = I, and X2 and ?/2 for i = 2. We 
shall assume that diagonal elements, even the ones obtained as the products of matrices, 
are local with respect to the coordinates xi and yi for the upper block and X2 and 1/2 
for the lower block, although they are in general bi-local. We will factorize a trace as 
Tr = Tr' tr2x2; and concentrate on a piece Tr'. There are three types of traces we have 
to consider, 
(i) A trace containing only diagonal elements, i.e., covariant derivatives: 



Tr [Dx^,Dy^]^ = <^o J dxidyi[Dr,^,D. 



2 



(7) 



where ctq = 5^ is a normalization factor and denotes the density of DO-branes on the 



(anti-)D2-brane worldvolume in the type IIA picture |T^ 



(ii) Two diagonal and two off-diagonal elements. An example is 

Tr'D^^YDy^X'' = a^ dxidyidx2dy2D^^Y{xi,yi;x2,y2)Dy^X*{xi,yi;x2,y2), (8) 

where * denotes the complex conjugation, and the same rule is applied to the other cases 

of this type. 

(iii) Four off-diagonal elements. An example is 

xX{xuyi; x's, y2)Y*{xi, yi, x'2, y'-^, (9) 

and similar results hold for the other cases. 

In order to see the appearance of tachyon modes in this system, it is instructive to 



compute the potential between M2 and M2 in the one-loop approximation. It can be 
done most easily in the background gauge with the gauge fixing term 

^o^ = -^^, J dt^Tr {d,Ao + t[B,, Y^]f . (10) 

Here the matrices F* denotes the fluctuations around the background, that is, X* = 
B^ _l_ Y^. Since we know that the tachyon modes would come from the open strings 
stretched between D2 and D2 in the corresponding type IIA picture, it suffices to consider 
the off-diagonal elements of the fluctuation matrices (and ghosts as well): 

Then the mass matrices for the off-diagonal fluctuations are given by 

Ml = -Ml = Mglosts = b' + H, (2 = 1, ■ ■ ■ , 7), (12) 

^($B,$9) = (^ _4^^ ^2 ^JJ j , (13) 

M^ = -7^6 + 78(Pi - P*) + 79(Qi + g*), (14) 

where H = {Pi — P^Y + {Qi + ^2)^' ^"^^ i^ describes the Hamiltonian of a charged 
particle moving on a two-dimensional plane in the presence of a uniform magnetic field 
perpendicular to the plane. Thus the energy level of the Hamiltonian H is identical to 
the Landau level, 4c(n + |) (n = 0, 1, ■ ■ ■). 



Now we can find that the tachyon develops when M2 and M2 approach to some critical 
interval, by diagonalizing the mass matrices for (8, 9)-directions: 

f/tM^ /6^ + 2c(2n + 3) \ n - — ('-' \ flS) 



The upper eigenvalue is the mass square of the matrix $ = -y^i^^ + «$^), and the lower 
is that of $ = 7f (*^^ —2$^). The mass square of $ becomes negative for the ground state 
when the interval b is shorter than a/2c, and thus the ground state of $ corresponds to a 



complex tachyon mode in the D2-D2 pair for b < \j1c |15, 16 



We can easily carry out the one-loop calculation, and the long range potential between 



Ml and M2 are reliably calculated in this approximation. The result is 

V{b) ~ -N r ^,-i/2,-.^. ;inh^c. _^_^ _ lr(5/2)^, (ig) 



/OS 2 sinh 2cs 4 b^ 

and the agreement with the supergravity result was found at least up to a numerical 
factor |l^, |16]. Speaking repeatedly, the potential ([T6| ) diverges for b < V2c, and it 
exhibits the tachyon instability of this system. A lesson from the one-loop computation is 
that the tachyon mode appears in the fluctuation $ when the distance between M2 and 



M2 becomes shorter than the critical value v 2c. 



3 The effective theory of the M2-M2 system 



We are now interested in the region b < v 2c. At the critical distance b = v 2c, a massless 
mode comes out and it might be related to the horizon of a black hole |T^ . We will mention 



this point later. For the region b < v 2c, a complex tachyon develops and thus it would 



annihilate the M2-M2 pair by its condensation. Moreover, as mentioned before, gravitons 
are supposed to be produced in this process as vortices on the membrane worldvolume. 



Let us discuss this mechanism in the case of the coincident M2-M2 pair, that is, 6 = 
case. For this purpose, it is enough to consider only the matrices X^ and X^, and to set 
the others to zero. The action in our concern is 



5'{8,9) 



^ fdt-Ti[X',X'^]\ (17) 



R^M^ 



The matrices X and X take the form 
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" 1, m^ P2 + ca^. }' \ i^y -Q2 - cay, ) ' ^^^> 

where a^;- and ay- denote the diagonal fluctuations. It is straightforward to show that the 
action (^) can be written, after setting $ = (deflned below the eq. (p!5D ), as 



%9) = -WmiI^^ 



— J dt (^Ic^ j dxidyidx2dy2\-\D^^'^ - F^y\^^\ 

f 1 ~ f 1 ~ 

+aoC^ J dxidyi-{F^y + F^yf + gqc'^ J dx2dy2-{F^y - F^yf 



CTn 



+ - 



2 



dxidyidx2dy2dx2dy2^{xi,yi; X2, y2)^*{xi,yi\ X2, 1/2) 

x$(xi, yi; X2, 2/2)'^*(2;i, 1/1; 4. 2/2) 
dx2dy2dxidyidx[dy[^* {xi, yi, X2, y2)^{xi,yi; X2, 2/2) 



X^*{x[,y[;x2,y2)^{x'i,y[;x2,y2) 



The covariant derivative D^ is defined by 
D, = 



Dy + iD^ 



Dy — dy +iAy +iay — [dy-^ — dy,^^ 



' 'f-{Axi — Ax2) + i[ctxi ~ Q-Z2J) 
h^yi + ^3/2) + HSi + S2)' 



and the field strengthes F^y and Fa;„ by 



-t xy 

Dx 



-i[Dx,D 



y\i 



xy 



-i[Dx,Dy], 



dx + iAx + idx = dx + i{Ax^ + A^J + i{ax^ + a 



X2J1 



D 



y 



Oy -\- lAy -\- l(Xy Oy + li^Ay-^^ ^1/2 J ' ^V*^?/! ^y2 ) ' 



(19) 

(20) 
(21) 
(22) 

(23) 
(24) 
(25) 



We shall restrict the configuration of a scalar field $ into the form, ^{xi,yi;x2,y2) 



^{x,y)J5{y/d^x)6{y/dQy), where x = |(xi— 0:2) and y = ^{yi + y2) are the center of mass 

(26) 



coordinates of the M2 and M2 system. Then the action (|19D reduces to 



S, 



4:an 



(8,9) 






dt / dxdy 
dt / dxdy 



— \D^\^-c^F m^+m^ + —F^+—F^ 



c 

4' ^-2^ 



c^^. 



xy 



\DM 



$1 



cH + c^ 



_c2f |<|.|2 I ^f2 I ^f2 
L- jxyl'i'l -r ^ Jxy ^ ^ Jxy 



C 

/^•ixy 



(27) 



where fxy = dxay - dyax{= Fxy - f ) and fxy = dxtty - dyax{= Fxy), and they are field 
strengthes for diagonal fiuctuations. When the tachyon field $ is at the lowest Landau 
level and the field strengthes fxy and fxy for diagonal fiuctuations are vanishing, the action 
(p7|) reduces to the sum of the tachyon potential | (|$p — c) — c^\ and the total amount 
c^ of tensions of M2 and M2. We can further rewrite the action (EgI) as 



S, 



4crn 



(8,9) 



R^M^ 



dt / dxdy 



2 ' ' I 2 ^ II 



—F: 



xy 



:28) 



and this provides us with a Bogomol'nyi-type equation p2 

-,2 



D.^ 



c 
~2 



' c\ ^y II ' 3^'j/ 



0. 



(29) 



The effective theory ( ^51) of the M2-M2 system is a U{2) gauge theory with the back- 
ground of a constant magnetic field, discussed by |^, § . 

At this point we would like to note that our membrane has a topology of two-torus 



p!3| . The periodicity of the (anti-)membrane is given by < Pj < y/ncN and < 



Qi < y/ncN {i = 1, 2) ||T5l. In terms of our specific representation (^) of Pi and Qi, the 



periodicity can be read off from the formulae Tyn/2In/2 = c"o / dxidyilN/2 and TrAr/2-^Af/2 = 
(Jo / dx2dy2lN/2, and the area At2 of torus, spanned by these coordinates, is equal to ncN, 
i.e., the same as that of {Pi, Qj)-planes. Let the periodicities of x and y he < x < Rx and 
< y < Ry respectively. Then R^. and Ry satisfy RxRy = A^^ = ncN. We would like to 
emphasize that the effective theory obtained in ( p8l) is not the Matrix theory compactified 
on T^. Our interpretation of the effective theory is similar to the one discussed in [|19| 



m 



which it is justified that the worldvolume theory obtained via the representation @ can 
be thought of as that on the membrane itself, not on the T-dualized torus. 

We would also like to remark that there is obviously a decoupled f/(l) gauge field 
in the effective theory. As pointed out in \W, TT], even the decoupled U{1) here should 



be eliminated in the M2-M2 annihilation. The author of |2^ proposed a resolution of 
this puzzle. It is not the same mechanism, but similarly we would like to point out that 
the decoupled f/(l) might be confined, that is, there is a mass gap and thus no massless 
vector in this U{1) theory, if we can justify that this U{1) is compact and the Polyakov's 
argument [^ goes through as well in the case of a two-torus. 

The background magnetic field gives a 't Hooft's twisted boundary condition p3|, |^. 
We shall take the following form of the background fields: 

A=(t -!) = ->- A=(t -X)=°- P°) 

Then the boundary conditions for these fields are given by 

A^{Rx,y) = -^nxd^,n-' + nxA^{o,y)n~\ 

Af,{x,Ry) = ^Qyd^Q-^ + QyAf,{x,0)Q-\ (31) 

where transition functions Qx and Qy are 

Qxix, y) = I2, Qyix, y) = ^^y^'^^l^ (32) 



and they satisfy a consistency condition given by 't Hooft | 



^Av = 0)fiy(x = Rx) = ^yix = 0)nxiy = Ry)i-ir. (33) 



Corresponding to these boundary conditions, we must have the boundary conditions for 
the other fields via 



X^=«'y(/2,, y) = fi.X^=«'y(0, y)n-\ X'=^'\x, Ry) 
and we can find that these conditions amount to 



fij,X*=^'9( 



X, 0)0-1 



y ' 



$(i?,,l/) = $(0,i/) 
l'(i?,,l/) = <l(0,i/) ' 



1 ^^y) 



<^{X,Ry 



^2iRyxlc^^^^ 0) 
g2iiJ,x/c^j^^ 0) 



(34) 



(35) 



4 Tachyon condensation and graviton production 

Now we shall show that there exist vortex-like solutions of the Bogomol'nyi-type equa- 
tion (pQ]). Let us express the 'Higgs' field $ by two real valued functions u and G as 
$ = exp[i(u + i6)]. From the first equation in (^) the f/(l) gauge field A^ \ a^ = 
\ [{Ar^ + a^) - i{Ay + ay)] is described by A^ + a^ = |9^(m + i6), where d^ = |(t?x - idy). 



log-i? 



' Z — Zl, 



log -(9 



U 27r 



Set the angular part 6 of the 'Higgs' field as 6 = z Y^k=i 

where the overline denotes the complex conjugation and B is a multivalued function.^ 

Then the second equation in (B^) reduces to 



c^d,d,u = -e" + 27Tc^ Y^^^i^' ^k) 



(36) 



k=l 



This is a Liouville-type equation discussed by Jackiw and Pi |25| , and general multi- vortex 
solutions are constructed at least on R^ ||2^ and an explicit one-vortex solution on T^ is 
given in [EH]. We, however, do not follow their constructions and adopt another strategy 



in order to handle the boundary conditions of vortex solutions manifestly. 

Now we divide the radial part u of the 'Higgs' field $ into a singular part Uq and an 
analytic part v. The singular function uq takes the form 



^o = J2{ lo^ 



fc=i 



^1 



Z — Zk 



2n 



1 2' 

[lm(2; - Zk)] 



2nT2 



(37) 



where T2 = Imr and this is exactly doubly periodic [^. The analytic function v obeys 
the equation 



c^dzdzV = c^ 



n 



AllTo 



-n 

fc=i 



^1 



Z- Zk 



271 



T 



e-2;^77[M^-^fe)]'g^-_ 



(38) 



^On R^, 8 may be chosen as —2 J22=i ^^si^ ~ ^k)- 



By the theorem 7.2 of Kazdan and Warner ||2^, one can show the existence of a solution 
for this equation, at least if -^^ is sufficiently small.Q 

Let us consider the boundary conditions for 'Higgs' field $. They can be easily read 
off from the angular part G of the 'Higgs' field: 

$(2 + 27rr) = e-*"'^+^™^""+*Sfc=i^"("-"*)<l>(z), (39) 

$(z + 27r) = e-^^"$(2). (40) 

Comparing these with (|35|) and noting that R^ = 2tt and Ry = 27it2 = cN/2 in this 
convention and that we can obviously replace x = Kez in ( P^ and ( P^f ) with -^ J2k=ii^ ~ 
Xk), one can find that those conditions coincide exactly by setting n = N{= even). 

At this point one may wonder if a smilar analysis goes through even for the case that 
$ 7^ and $ = 0, that is, regardless of the existence of the tachyon.Q It is almost true. 



The Bogomol'nyi-type equation (|29D becomes -Dg$ = where D^ = —iD^ + Dy, and 



Y^xy + I^P = 0, and thus we can solve it similarly. There is, however, a significant 
difference. In this case the angular part ffips its sign and we obtain the boundary 
conditions for $ with wrong sign. Thus we cannot construct vortex solutions for $ with 
the correct boundary conditions. 

Now let us discuss the tachyon condensation in this system. The effective theory is 
a pure U{2) gauge theory and thus it seems that the tachyon cannot condense. We can, 
however, show that it eventually does. Integrating (|38|) over T^, we can rewrite it as 

c^ At2 = / dzdmi"^. (41) 

From the above argument, we have N/Attt2 = 1/c. Therefore we can conclude that the 
average value of |$| is equal to ^/c. It is anticipated by the form of the tachyon potential 
{(|$p — c)^ — c^} which we observed before, and this shows the tachyon condensation in 



the M2-M2 system. 

Now the magnetic field F^y on the worldvolume of the coincident M2-M2 pair is given 
by F^y = J^ 2dzdzV, and thus the magnetic fiux is 



^y 2ttt2 



I dzdzFxy = 2'kN (= 27r • 2ao ■ At2). (42) 

^ 1 



^Their theorem is the foUowing. Consider a differential equation dzdgv = C — He^, on Ricmann 
surfaces, where C is a positive constant and H is an analytic function. Then there exists an analytic 
solution V, if and only if H is positive somewhere and C < C+{H), an upper bound depending on H . But 
in our case we can expect that n/47rT2 does not need to be small, since generic solutions of the Liouville 
equation ( p6[) is known and it ensures very likely the existence of a solution of eq. ( pq) . 

^Throughout this paper we are assuming c > 0. 



This corresponds to the A^- vortex solution on T^, in which each vortex is sitting at one 
of the zeros Zk of the 'Higgs' field. As is well-known, a magnetic flux on the membrane 
worldvolume couples to a RR one-form via a Chern-Simons coupling and thus gives a 
Z^O-brane charge in the corresponding type IIA picture. The number of DO-branes corre- 
sponding to the above flux is equal to A^, since the value 2'kN of the flux is just twice the 
one {J dxidyi- = J dx2dy2- = nN) on each ( ant i-) membrane which is constituted from 
N/2 DO-branes. Therefore the process we obtained will be 



M2 (+ N/2 M-mom) + M2 (+ N/2 M-mom) ' '^ ) gravitons (+ A^ M-mom), (43) 

where "M-mom" indicates the M-momentum, and this process in fact conserves the M- 
momentum. We will justify below that vortices can be interpreted as gravitons, as denoted 
in the r.h.s. of the above equation. 

Now let us show that vortices can be thought of as gravitons in Matrix theory. The 



bound for the energy is exactly zero |^. This indicates that vortices in this system are 
massless particle, and thus they can be regarded as gravitons in M-theory. It, however, 
seems somewhat strange that a massless particle does not carry any energy. But remember 
that the matrix theory Hamiltonian is of the form H = (P^ + M^)/2P_. Hence in matrix 
theory the Hamiltonian of a massless particle without transverse momentum is exactly 
zero, which accords with the above result. We would like to note that anti-gravitons 
cannot appear in this system, since as we have shown the magnetic flux is positive. Again 
it is in accordance with the interpretation of Matrix theory |13|. A graviton of this 
type carries the M-momentum in the form of a magnetic flux / dxdyFxy, as mentioned 
above, which may be reminiscent of the membrane scattering with M-momentum transfer 
3^|3l[|. A vortex solution is, however, a static solution, as is different from an instanton 



solution which can provide scattering processes with M-momentum transfer. Hence a 
vortex solution cannot describe a transfer of M-momentum. 

We would also like to remark that, at the minimum of the tachyon potential, the sum of 



the tension of M2 and M2 is exactly cancelled by the energy density of the tachyon field, 
as discussed in 1^,0. We can observe it from the action (P7|): (here in unit ao/R^M^ = 1) 

Em2 + Ej^= -^Ti[B',B'f = (^Jdxidyi + J dx2dy2^ ^c' = AJdxdyc\ (44) 

^tachyon = A j dxdy [m'' - cf - c^] ^^ - A j dxdy c\ (45) 

This indicates the restoration of supersymmetry ||^ . Now let us discuss the supersymmetry 
of this system. The initial state, the M2-M2 pair, breaks all of the supersymmetry. As we 
showed above, it is annihilated via the tachyon condensation and gravitons are produced 

9 



as vortices on the membrane worldvolume. The final state, gravitons, is a massless state, 
and thus it must preserve a half of the supersymmetry. Let us look at the fermion 
transformation law under the supersymmetry transformations: 

5^ = [OtX^^i - '-[X\ X^]7,,) e + e'. (46) 



For a vortex solution we can find that 

2 I-' xj/ T^ -f xy) \^ 
c 

V2 



6^ = ' 



2\ M^^^y t\^{F,y-F,y)-m 



e . 




This shows that a vortex solution indeed preserves a half of the supersymmetry, as we 
expected. Note also that gravitons in Matrix theory are described by a background of 
commuting matrices. The above result (^) is nothing but the statement that X^ and 
X^ are commuting for a vortex solution. Again this agrees with the interpretation of a 
vortex as a graviton in Matrix theory. 

5 Discussions 

Now we would like to discuss future directions of this work. The authors in |jl8| pointed 
out the similarity of the tachyon instability to the absorption of matters by a non-extremal 
black hole. They argued that this similarity leads to the identification of the horizon of 
a non-extremal black hole with a critical point (corresponding to a point b = \/2c in 



the M2-M2 system) where a massless mode comes out. As the horizon has a light-cone 
structure, the appearance of a massless mode matches this feature. Along the line of 



their arguments, it would be interesting to pursue the following possibility in the M2-M2 
system. We will indicate the corresponding picture in terms of a black hole in parenthesis: 



Set M2 and M2 infinitely separated, that is, b(t = — oo) = oo. The attractive force 



between M2 and M2 makes them approach up to a critical point (the horizon), b(t) = v 2c, 
where a massless mode appears, and subsequently they across this point (the horizon) 
to exhibit the tachyon instability. Then there appear vortices (gravitons radiated by a 
black hole) which are time-dependent through b{t) and thus carry transverse momenta in 



the region, < b{t) < v2c- We expect that M2 and M2 gradually collapses to a point, 
b{t = oo) = 0, and there remain only vortices (gravitons) as discussed above. In this 



senario the initial energy of M2 and M2 seems to disappear, for a vortex at the final state 
does not carry any light-cone energy as mentioned before. But Matrix theory is unitary 



and thus the total energy of M2 and M2 must be radiated (the Hawking radiation) by 
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a certain mechanism in which vortices (gravitons), in the region < b{t) < v2c, carry 
away the initial energy. 

Another direction is to study stable non-BPS states [§, |, |^, |, || in Matrix theory. To 
do this, we need certain orbifold or orientifold projections to make a complex tachyon real 
so that there can exist a kink solution on the brane - anti-brane pair. As an example it 
is straightforward to apply our analysis to a Matrix theory on an orientifold background 



1 3^ . We will report some progress in this direction elsewhere ||3^. It also seems possible 



to discuss stable non-BPS states in type I theory p, by making use of the USp{2k) 



matrix model ^A 



Finally we would like to mention that our analysis is applicable to more general sys- 



tems, such as general brane - anti-brane pairs in the IIB matrix model [^ as well as 
those in Matrix theory. 
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